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ABSTRACT 

We introduce a new notion of "matrix potential" to nonlinear optical systems. In 
terms of a matrix potential g, we present a gauge field theoretic formulation of the 
Maxwell-Bloch equation that provides a semiclassical description of the propagation 
of optical pulses through resonant multi-level media. We show that the Bloch part of 
the equation can solved identically through g and the remaining Maxwell equation 
becomes a second order differential equation with reduced set of variables due to the 
gauge invariance of the system. Our formulation clarifies the (nonabelian) symmetry 
structure of the Maxwell-Bloch equations for various multi-level media in association 
with symmetric spaces G/H. In particular, we associate nondegenerate two-level 
system for self-induced transparency with G/H = SU{2)/U{1) and three-level A- 
or l^-systems with G/H = SU{3)/U{2). We give a detailed analysis for the two- 
level case in the matrix potential formalism, and address various new properties of 
the system including soliton numbers, effective potential energy, gauge and discrete 
symmetries, modified pulse area, conserved topological and nontopological charges. 
The nontopological charge measures the amount of self-detuning of each pulse. Its 
conservation law leads to a new type of pulse stability analysis which explains nicely 
earlier numerical results. 
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1 Introduction 

Since the invention of the laser, much progress has been made in understanding 
nonlinear interactions of radiation with matter which made nonlinear optics a fast 
developing and independent field of science. Recently, the interaction of laser lights 
with a multi-level optical medium has attracted more attention in the context of 
lasing without inversion and electromagnetically induced transparency (EIT) 
101 . Laser light in general is expressed in terms of a macroscopic, classical elec- 
tric field which interacts with microscopic, quantum mechanical matter. Unlike 
classical electrodynamics, the electric scalar potential and the magnetic vector po- 
tential do not appear to replace electromagnetic fields in nonlinear optics. Instead, 
the electric field itself, with appropriate restrictions to accomodate specific physi- 
cal problems, plays the role of a fundamental variable which renders the problem 
lacking a field theoretic Lagrangian formulation. Of course, one could setup the 
problem in the most general QED Lagrangian framework with the conventional 
potential variable A^, but the nonlinearity of interactions and various approxima- 
tion schemes involved make the use of potential meaningless. For instance, the 
Maxwell-Bloch equation which governs the interaction between radiation and mat- 
ter takes a nonrelativistic, semiclassical limits of QED together with slowly varying 
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envelope aprroximation (SVEA) and/or rotating wave approximation (RWA). Vari- 
ables of the Maxwell-Blocli equation are given by the envelope functions of electric 
fields, and the components of the density matrix or the probility amplitudes for each 
atomic level occupations. Thus, all previous works have focused on the study of the 
Maxwell-Bloch equation itself, without making any reference to the Lagrangian and 
potential variables. However, there exists one notable exception. In the case of 
nondegenerate two-level atoms, McCall and Hahn |Q have shown that lossess prop- 
agation of light pulses, the phenomenon of self-induced transparency (SIT), can be 
explained in terms of a potential-like variable 6{x), the time area of a suitably cho- 
sen electric field, which obeys the area theorem. Under certain circumstances, the 
system can be described by an effective potential variable <f{x,t) which satisfies the 
well-known sine-Gordon equation. In this case, the 1-soliton of the sine-Gordon the- 
ory is identified with the 27r-pulse of McCall and Hahn. The cosine potential term 
becomes proportional to the microscopic atomic energy, and the stability of the 27r- 
pulse is explained through the topological charge conservation law. Recently, the 
quantum sine-Gordon theory has been also applied to the Maxwell-Bloch equation 
and quantum optics with interesting results 0. However, one serious drawback of 
the sine-Gordon approach to the Maxwell-Bloch system is its oversimplification. In 
the sine-Gordon limit, frequency detuning and frequency modulation effects are all 
ignored and microscopic atomic motions (inhomogeneous broadening) are not taken 
into account. Also, the model is limited only to the nondegenerate two-level case 
while many recent interesting applications are based on the multi-level (three-level 
and higher) and possibly degenerate systems. In an earlier work 0], we have shown 
that even the nondegenerate two-level system should be described by the complex 
sine-Gordon equation. This generalizes the sine-Gordon equation by including a 
phase degree of freedom which accounts for frequency modulation effects. We have 
also shown that a more general framework can be given by a 2 x 2 matrix potential 
g and its Lagrangian formulation. This allowed us to incoporate frequency detun- 
ing and external magnetic fields. Until now, the sine-Gordon theory was the only 
available field theory for the Maxwell-Bloch system and therefore all analytic works 
beyond the simplest two-level case have resorted to the Maxwell-Bloch equation, 
finding soliton type solutions through the inverse scattering method in integrable 
cases (for a review, see ^ and other references therein). Following the pioneering 
work of Lamb P], Ablowitz, Kaup and Newell have extended the inverse scattering 
formalism to include inhomogeneous broadening and obtained exact solutions 
In accordance with the area theorem, these solutions show that an arbitrary initial 
pulse with sufficient strength decomposes into a finite number of 2tt pulses and Ovr 
pulses, plus radiation which decays exponentially. Extensions to the degenerate as 
well as the multi-level cases have been also found resulting more complicated soliton 
solutions 0, |T0|, [n], g. 



In this paper, we introduce a new matrix potential variable g to nonlinear opti- 
cal systems described by (integrable) Maxwell-Bloch equations, and present a com- 
pletely different type of analysis of the Maxwell-Bloch equation based a field theory 
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formulation through g. We show that the Bloch part of the equation can be solved 
identically in terms of g and the remaining Maxwell part becomes a second order 
differential equation in g. This is compared with the linear case of electromag- 
netism where the curl-free condition is solved in terms of a scalar potential <f and 
the remaining Gauss equation changes into the second order differential equation in 
if. The field theory action for the second order Maxwell equation in g is provided 
by a sigma model-type action which combines the so-called "the 1+1- dimensional 
G / H-gauged Wess-Zumino-Novikov-Witten action" with an appropriately chosen 
potential energy term. This work which generalizes the earlier work on the two- 
level case to the multi-level cases uncovers many new features of the problem. 
In particular, our formulation clarifies the hidden (nonabelian) group structure of 
the multi-level Maxwell-Bloch equation in association with symmetric spaces G/H. 
For instance, nondegenerate two-level system of self-induced transparency is asso- 
ciated with G/H = SU{2) /U{1) while three-level A- or l^-systems are associated 
with G/H = SU{3)/U{2). These nonabelian group structures are shown to arise 
from the probability conservation law of a density matrix and also from the selection 
rules in relevant dipole transitions. In general, the number of degrees of freedom 
for the Maxwell equation (those of electric field components) is smaller than that 
of the matrix potential g belonging to the group G. We show that these residual 
degrees can be removed by imposing constraints on g through "gauging" the action 
so that the action possesses the if-vector gauge invariance. The gauge transforma- 
tion, however, is shown to receive physical meaning at the atomic level. That is, it 
accounts for the effects of frequency detuning and external magnetic fields. We show 
that inhomogeneous broadening can be also incoporated into the matrix potential 
formalism. 

In order to demonstrate the power of our matrix potential approach, we make 
a detailed analysis of optical pulses. This shows that the matrix potential not only 
leads to a deeper understanding of optical pulses, but it also provides new solutions, 
new conserved charges and symmetries. In particular, a new type of stability anal- 
ysis is made which generalizes the area theorem to a certain extent. Specifically, we 
clarify the topological nature of solitary pulses through the effective potential energy 
term and its degenerate vacua. We define the topological soliton number accord- 
ing to the group structure of the system and show that a sohtary pulse for certain 
multi-level cases, e.g. the degenerate three-level case, carry more than one soliton 
numbers. Also, we show that 27r pulses can be nontopological carrying a nontopo- 
logical charge. A nontopological soliton is interpreted as a "self-detuned" 27i pulse 
and the nontopological charge is shown to measure the amount of frequency self- 
detuning. The conservation laws of the topological and the nontopological charges 
are shown to prove the stability of pulses. In particular, we prove the stability of 27r 
pulses against small fluctuations. This explains nicely the frequency pulling effect in 
the presence of frequency detuning which has been predicted earlier by a numerical 
work. 

Our matrix potential formalism also allows a systematic understanding of various 
symmetry structures of the Maxwell-Bloch equation. We show that infinitely many 
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conserved local integrals resulting from the integrability of the equation can be 
obtained in a general group theoretic framework of symmetric space G/H. This 
enlarges previously known results in the case of the two-level system and provides 
new conserved charges in other multi-level cases. More importantly, our field theory 
reveals new types of symmetries; i) global gauge symmetry, ii) global [/(l)-axial 
vector symmetry, iii) chiral symmetry and iv) dual symmetry. We show that global 
gauge symmetry can be used to generate simulton solutions systematically. Global 
f/(l)-axial vector symmetry gives rise to the nontopological charge via the Nocthcr 
method. Chiral and dual symmetries are discrete symmetries and they generate 
new solutions from a known one. In particular, dual symmetry relates the "bright" 
soliton with the "dark" soliton of SIT. Finally, we show that the matrix potential is 
useful in understanding the inverse scattering method itself. The potential variable 
g reveals the group structure of the inverse scattering method and we construct 
explicitly soliton solutions for various cases. 

The plan of the paper is the following; in Sec. 2, we present a field theory 
formulation of the Maxwell-Bloch equation. The area theorem and the sine-Gordon 
field theory limit are briefiy reviewed and an extension to the complex sine-Gordon 
field theory is made in Sec. 2.1. In Sec. 2.2, a matrix potential formalism is presented 
and a general action principle is found for the Maxwell-Bloch equation for arbitrary 
multi-level systems. In Sec. 2.3, inhomogeneous broadening is also incoporated into 
the matrix potential formalism. Section 3 deals with explicit examples of various 
multi-level systems. Specific group structures and gauge fixing for each systems are 
identified. In Sec. 4, we explain new features of optical pulses in our matrix potential 
formalism. In Sec. 4.1, topological properties of pulses are analyzed through the 
effective potential energy and its degenerate vacua and also by defining topological 
soliton numbers. In Sec. 4.2, nontopological solitons are introduced and interpreted 
as self-detuned pulses. In Sec. 4.3, a new analysis of pulse stability is made in 
terms of newly found nontopological charges. Section 5 deals with symmetries of 
the system. Infinitely many conserved charges are constructed systematically for the 
general multi-level systems in Sec. 5.1. Global gauge symmetries are explained in 
Sec. 5.2 and the chiral and the dual symmetries are explained in Sec. 5.3. Finally, 
Sec. 6 is a discussion. 

2 Field theory for the Maxwell-Bloch equation 

The multi-mode optical pulses propagating in a resonant medium along the x-axis 
are described by the electric field of the form, 

m 

E = ^ Ei{x, t) exp i{kix - wit) + c.c. (2.1) 
1=1 

where ki and wi denote the wave number and the frequency of each mode and the 
amplitude vector E; is in general a complex vector function. The governing equation 
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of propagation is the Maxwell equation, 



::t7si^)e = -t^ M^trpd. (2.2) 



(9x2 ^2 ^^2 ^2 ^^2 

On the right hand side, electric dipole transitions are treated semiclassically. d is 
the atom's dipole moment operator and the density matrix p satisfies the quantum- 
mechanical optical Bloch equation 

in{--rV—)p=[{H,-'E-d), p]. (2.3) 

Hq denotes the Hamiltonian of a free atom and v is the x-component of the velocity 
of the atoms. In general, we make a slowly varying envelope approximation (SVEA) 
for the Maxwell-Bloch system where the amplitudes vary slowly compared to 
the space and time scales determined by ki and wi. Under SVEA, the Maxwell- 
Bloch equation becomes a set of coupled first order partial differential equations 
for the amplitudes E^ and the components of the density matrix. Explicit expres- 
sions of the Maxwell-Bloch equation for several multi-level cases are given in Sec. 
3. Thus, the Maxwell-Bloch equation provides an effective, semiclassical description 
of light-matter interaction using the amplitudes E^ as dynamical variables. Unlike 
the linear case, it is quite difficult, if not impossible, to introduce a potential vari- 
able instead of E; due to the nonlinearity of the interaction and the approximation 
involved. Lacking a potential variable causes the physical system to be described 
only by the equation of motion, not by an action principle. Consequently, a field 
theoretic formulation is lacking in the problem of pulse propagation. However, when 
pulses propagate in a resonant, nondegenerate two- level atomic medium with inho- 
mogeneous broadening, McCall and Hahn have introduced an effective potential-like 
variable, and in terms of which shown that an arbitrary pulse evolves into a coherent 
mode of lossless pulses [Q. This phenomenon, known as self-induced transparency, is 
also observed in more general, degenerate and/or multi- level atomic media. Specif- 
ically, McCall and Hahn have shown that when the dimensionless pulse envelope 
function E is assumed to be real and the time area of 2E, 

e{x) = 2 dt E, (2.4) 



is an integer multiple of 2n ( 2mT pulse), then the pulse propagates without loss 
of energy. Otherwise, due to inhomogeneous broadening the pulse quickly reshapes 
into a 2n7i pulse according to the area theorem, 

^ = -asin0(x) (2.5) 
dx 

for some constant a. The proof of the area theorem can be done by making use 
of inhomogeneous broadening and the Maxwell-Bloch equation. In the absence of 
inhomogeneous broadening, the system was shown to admit the sine- Gordon field 
theory formulation. 
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2.1 The sine- Gordon limit 

The Maxwell-Bloch equation for the nondegenerate two-level case can be written in 
a dimensionless form by 

dE + 2p<P> = 
dD - E*P - BP* = 

dP + 2i^P + 2ED = (2.6) 

where /3 is a coupling constant and ^ = w — wo,d = d/dz, d = d/dz , z = 
t — x/c, z = x/c. The angular bracket signifies an average over the spectrum 
distribution /(^) as given by 

/oo 
(■■■)fim- (2.7) 
-oo 

The dimensionless quantities E, P and D correspond to the electric field, the polar- 
ization and the population inversion through the relation, 

E = -iE-etod/Ven 

P = -puexp[-t{kx-ujt)]/AktoNof{0 

D = -(p22-pii)/8HoiVo/(0 (2.8) 

where e specifies the linear polarization direction, to is a time constant and A^o is 
related to the stationary populations of the levels.^ In order to understand the 
structure of 2mT pulses better, we impose further restrictions such that the system 
is on resonance = 0), frequency modulation is ignored {E being real) and inho- 
mogeneous broadening is absent (/(O = ^iO)- Under such restrictions, we could 
introduce an area function ip{x,t) defined by 

/■* 

(p{x,t) = Edt , (2.9) 

J —oo 

which, in the limit t — > oo, agrees with 6{x)/2 in ( p.4|) . In terms of ip, the SIT 
equation reduces to the well-known sine-Gordon equation, 

ddifi -2psm2ip = 0, (2.10) 

when we make consistent identifications; 

E = E* = dip , < P >= P = - sin2ip , < D >= D = coii2ip . (2.11) 

This sine-Gordon equation arises from the action 

S = ^ f {d^d^ - 2/3 cos 2y?) . (2.12) 
27r J 



^For the details of constants, we refer the reader to Ref. 
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The periodic cosine potential term exhibits infinitely many degenerate vacua. It 
gives rise to soliton solutions which interpolate between two different vacua. This 
shows that the 2mT pulse can be identified with the topological n-soliton solution of 
the sine-Gordon equation. The electric field amplitude E, now identified with d(f, 
receives an interpretation as a topological current. Note that the area function (p 
is different from the conventional scalar or vector potentials of electromagnetism. 
Nevertheless, it is remarkable that the potential energy cos2(/? of the sine-Gordon 
Lagrangian can be identified with the population inversion D which represents the 
atomic energy. Also the Lorentz invariance, which was broken by SVEA, re-emerges 
in the sine-Gordon field theory after the redefinition of coordinates. The identifica- 
tion of the atomic energy with the cosine potential term shows that 2n7r pulses are 
stable against finite energy fiuctuations due to the conservation of the topological 
number n. 



Though the sine-Gordon theory provides a nice field theory for the nondegenerate 
two-level system, it is too restrictive for real applications. The presense of frequency 
modulation in pulses, for example, require that E should be complex. Therefore, 
in this case E can not be simply replaced by a real scalar field and the sine- 
Gordon limit is no longer valid. Also, inclusion of frequency modulation invalidates 
the area theorem. However, through the inverse scattering method, it has been 
found that solitons do exist even in the case of complex E [§]. This suggests that 
a more general field theory of SIT than the sine-Gordon theory could exist which 
takes care of a complex E. Recently, we have shown that this is indeed true and 
the field theory which includes both the frequency detuning and the modulation 
effects is the so-called "complex sine-Gordon theory" ||^. This generalizes the sine- 
Gordon theory as follows; assume that E is complex and the frequency distribution 
function of inhomogeneous broadening is sharply peaked at ^, i.e. /(^ ) = 6{C, — ^) 
for some constant ^. Introduce parametrizations of E, P and D, which generalize 
parametrizations in Eq. ( |2.11| ), in terms of three scalar fields p, 6 and t], 



E = e^(e-2r?)(29r/^^ - idcp) , P = ^e'^'^-^r,) ^■^2ip , D = cos2ip . (2.13) 
sin (f 

These parametrizations consistently reduce the two-level Maxwell-Bloch equation 
(|2.6| ) into a couple of second order nonlinear differential equations known as the 
complex sine-Gordon equation; 

— cos — 

dd(p + A — :T^dr]d7] -2psm2ip = 
sm ip 

2 - 

ddrj {drjdp + drjdip) = (2-14) 

sin 2(p 

and a couple of first order constraint equations, 

2 cos^ ifdr] — sin^ pdO — 2^ sin^ = 

2 cos Vt^^ + sin V^^^ = 0. (2.15) 



8 



Note that the complex sine-Gordon equation reduces to the sine-Gordon equation 
when frequency modulation is ignored so that 77 = 0, Q = -njl and the system is 
on resonance = 0). This reduction is consistent with the original equation since 
solutions of the sine-Gordon equation consists a subspace of the whole solution space. 
The complex sine-Gordon equation was first introduced by Lund and Regge in 1976 
in order to describe the motion of relativistic vortices in a superfluid |13[, and also 



independently by Pohlmeyer in a reduction problem of 0(4) nonlinear sigma model 
. This equation is known to be integrable and soliton solutions generalizing those 



of the sine-Gordon equation have been found. These issues will be considered in later 
sections in a more general context. The Lagrangian for the complex sine-Gordon 
equation in terms of and ri is given by 

S = ^ f d<pd<p + 4 cot^ <^dr]dr] - 2/5 cos 2v7. (2.16) 

This Lagrangian, however, is singular at (f = rnr for integer n which causes difficul- 
ties in quantizing the theory. Also, besides the complex sine-Gordon equation, the 
two- level Maxwell-Bloch equation comprises the constraint equation ( p. 151 ). Thus 



the Lagrangian (|2.16|) does not quite serve for a field theory action of the two- 
level system. In fact, the singular behavior of the Lagrangian ( p.l6[ ) is an artifact 



of neglecting the constraint equation. This, as well as the rationale of the above 
parametrizations, can be seen most clearly if we reformulate the Lagrangian to in- 
clude the constraint in the context of a matrix potential and a gauged nonlinear 
sigma model as explained in the next section. 

2.2 Matrix potential formalism 

In order to construct a field theory action of the Maxwell-Bloch equation in terms 
of potential variables and also find a way to extend to more general multi-level 
and degenerate cases, we first note that the optical Bloch equation admits an in- 
terpretation of a spinning top equation as in the case of the corresponding mag- 
netic resonance equations |]T3|]. Denote real and imaginary parts of E and P by 
E = En + iEi, P = Pr + iPj. Then, the Bloch equation in Eq. (|2.6| ) can be 
expressed as 

dS = nxS (2.17) 

where 5* = (Pr, P/, D), VL= {2Ej, — 2Er, — 2^). This describes a spinning top 
where the electric dipole "pseudospin" vector S precesses about the "torque" vector 
fl. This clearly shows that the length of the vector S is preserved, 

\S\'^ = P^ + Pf + = 1. (2.18) 

The length equals unity due to the conservation of probability. The remaining 
Maxwell equation in Eq. ( p.6|) determines the strength of the torque vector. If 
Pj = 0, we may solve ( |2.18D by taking Pji = — sin 2ip and D = cos 2ip and also ( p.l7| ) 
by taking E = dip as given in ( p.ll|) . Then, the Maxwell equation becomes the sine- 



Gordon equation as before. This picture agrees with the conventional interpretation 
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of the sine-Gordon theory as describing a system of an infinite chain of pendulums. 
In order to generalize the sine-Gordon hmit to the complex E and P case, we make 
a crucial observation that the constraint in Eq. (|2.18|) can be solved in general in 
terms of an SU{2) matrix potential variable g by 



D P 

P* -D 



-'a,g, (2.19) 



where as = diag(l, —1) is the Pauli spin matrix. By taking the determinent, one can 
check that Eq. ( p.l8|) is automatically satisfied. Also, note that g'^cTsg is invariant 
under the "chiral f/(l)-transformation" , 

9 ^ e'^'^'g (2.20) 

for any function /. Thus, g^^cr^g parameterizes SU{2)/U{1) instead of SU{2). 
Since Eq. ( [^.18| ) is automatically solved in this SU{2) parametrization, the number 
of independent variables, which is dim(5'f/(2)) — dim(?7(l)) = 3 — 1 = 2, agrees with 
that of the vector {S) parametrization. Moreover, we have an identity 

dig-'cx,g) = [g-'a,g, g-'dg], (2.21) 

where the bracket denotes a commutator. Note that this becomes precisely the 
Bloch equation if we make an identification, 

9-'dg + R=^'j^^ :|), (2.22) 

where R is an anti-Hermitian matrix commuting with g^^cr^g which will be deter- 
mined later. Thus, we have solved the Bloch equation through the matrix potential 
g up to the identification in Eq. (|2.22|) . The identification is consistent since both 



sides are anti-Hermitian matrices. The off-diagonal part of the r.h.s. is simply re- 
naming the component variable by E whereas the constant diagonal part imposes 
a constraint on the variable g. This constraint, however, can be satisfied by an 
appropriate chiral [/(l)-transformation in Eq. (|2.20|) . Thus, the matrix potential g 
is made of two independent variables and one variable satisfying the constraint. In 
the following, we show that the Maxwell equation can be expressed in terms of two 
independent variables only, decoupling completely from the constraint variable. In 
this regard, our matrix potential g resembles the scalar potential (p in electrostat- 
ics where if solves the curl-free condition, V x = 0, identically and changes the 
Gauss equation into the Poisson equation. In our case, the Schrodinger equation 
plays the role of the curl-free condition and the Maxwell equation, the counterpart 
of the Gauss equation, changes into a second order nonlinear differential equation. 
In order to see this, observe that the Maxwell equation can be expressed also in 
terms of g only, 

d{g-'dg + R) = 'f^) = (3[a,, g~'a,g]. (2.23) 
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Thus, we have successfully expressed the SIT equation in terms of the potential 
variable g up to an undetermined quantity R. As we will show, R is determined 
by requiring an action principle for the Maxwell equation in terms of g. Since g 
is constrained, we need a Lagrange multiplier for the constraint. In order to help 
understanding, we assume for a moment that i? = and the system is on resonance 
(i^ = 0). Then, the equation of motion (|2.23|) arises from a variation of the action 



•S* - SwzNwia) - 'S'pot + 'S'const (2.24) 
with the following variational behaviors; 

8gSwzNW = ^ J dzdzTT{d{g'^dg)g~^Sg) 
P f 

^gSpot = ^ y dzdzTi^as, g'^asg ]g'^6g). (2.25) 

The action SwzNw{g) is the well-known SU{2) Wess-Zumino-Novikov-Witten func- 
tional, 

SwzNwig) = f dzdzTT (g-^dgg'^dg) - f Tr {g'^dg A g'^dg A g'^dg) , 

Ait jt. IZtt Jb 

(2.26) 

where the second term on the r.h.s., known as the Wess-Zumino term, is defined on 
a three-dimensional manifold B with boundary S = dB and g is an extension of a 
map (yf : S — > SU{2) to B with = g ljTB[- The potential term can be easily 
obtained by 

5'pot = ^ J dzdzTrigasg-'^as). (2.27) 

Finally, the constraint requires vanishing of the diagonal part of the matrix g~^dg 
which can be imposed by adding a Lagrange multiplier term •S'const action 

-^const = 7^ / dzd-zTi{\a^g-^dg). (2.28) 

The Lagrange multiplier A, however, induces a new term to the equation of motion 
by 

^gSconsi = ^ I dzd-zTi{{-d\a^ + [Aag, g-^dg])g-^5g), (2.29) 

which seems to spoil our construction of a field theory. This problem can be resolved 
beautifully if we introduce a gauge symmetry and make the action ( |2.24| ) to be 
"vector gauge invariant" . This can be done by replacing the constraint term with a 
"gauging" part of the Wess-Zumino-Novikov-Witten action, 

5" = Swzw{g) — 'S'pot + 'S'gauge (2.30) 
^-gauge = ^ / M-Adgg-' + Ag~'dg + AgAg-' - AA) (2.31) 
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where the connection fields A, A gauge the anomaly free subgroup U{1) of SU{2) 
generated by the Pauli matrix a^. They introduce a ?7(l)-vector gauge invariance 
of the action where the f/(l)-vector gauge transformation is defined by0 

g^h-^gh , A^ h-^Ah + h-^dh , A ^ h'^ Ah + h-^dh (2.32) 

where h = exp{ifas) for some scalar function /. Owing to the absence of kinetic 
terms, A, A act as Lagrange multipliers which result in the constraint equations 
when A and A are integrated out. The action in Eq. (|2.24|) may be understood as a 
gauge fixed action with the choice of gauge where A = 0, A = Xa^. The main reason 
for introducing a gauge transformation and a gauge invariant action is twofold. It 
first shows that the equation of motion resulting from the variation 6gS = 0, 

d{g-'dg + g-'Ag) + [A, g-'dg + g'^Ag] - dA = P[as, g-'a^g], (2.33) 

is also gauge invariant and gives rise to the gauge invariant expression of the Maxwell 
equation. Comparing Eq. (|2.33| ) with Eq. (|2.23| ), we see that R = g~^Ag and the 



relevant gauge choice is A = i^a^, A = due to the constraint in Eq. ( p.43| ) .p| 
The [/(l)-vector gauge invariance of the Maxwell equation implies that the Maxwell 
equation is independent of specific gauge choices. Thus, it decouples from the U{1) 
scalar field which saturates the constraint condition and becomes a couple of second 
order nonlinear differential equations in two local variables. In the next section, 
this is shown clearly by an explicit parametrization of g and the resulting Maxwell 
equation is shown to be equivalent to the complex sine-Gordon equation given in Eq. 
2.14). The second reason is that gauge transformation incoporates beautifully the 



frequency detuning effect through specific gauge fixing. In Sec. 3, external magnetic 
fields are also incoporated through gauge fixing. 

Our field theory for the Maxwell-Bloch equation is not restricted to the two-level 
case. In fact, the group theoretical formulation through the matrix potential g allows 
an immediate extension to the multi-level cases. We may simply replace the pair 
SU{2) D t/(l) by G D if for any Lie groups G and H and obtain the G/if-gauged 
Wess-Zumino-Novikov-Witten action {Swznw + 5'gauge) where A and A gauge the 
subgroup H of For a general pair of G and if, the expression for the potential 
which preserves integrability can be given by ||19, 20 1, 



^pot = ^ / "^^{gTg-'T) (2.34) 



4 



Note that here we are using the vector [/(l)-transformation instead of the chiral one 



as in Eq. (2.20). This causes the matrix g~^(J^g to transform covariantly under the gauge 
transformation. 

^One can always choose such a gauge due to the flatness of A and A as in Eq. (2.44). 



®This action is known to possess conformal symmetry and has been used for the general 
G/iJ-coset conformal field theories The potential energy term (2.27) breaks conformal 



symmetry. Nevertheless, it preserves the integrability of the model given by (2.30) where 
G/H = SU{2)/U{1), and this model has been used in describing integrable perturbation 
of parafermionic coset conformal field theories [pSl, E9| . 
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where T and T are constant matrices which commute with the subgroup H, i.e. 
[T, /i] = [T, }i\ = 0, for h & H. This makes the potential term vector gauge invari- 
ant .[| As we will see later, physically interesting cases all correspond to a special 
type of symmetric spaces G/H, known as Hermitian symmetric spaces ||2^, where 
the adjoint action of T defines a complex structure on G/H. 

Now, we define the field theory action for the Maxwell-Bloch equation by 

Smb = SwzNw{g)+^ I -TTi-Adgg-^+Ag-^dg+AgAg-^-AA)-^ f 1i{gTg'^T). 

ATI J zn J 

(2.36) 

This action is of course restricted to the integrable cases which require specific fine 
tuning of coupling constants. However, the concept of matrix potential g is valid 
irrespective of the integrability of the model and the field theory formulation can 
be extended to more general, nonintegrable cases too. In this paper, we will restrict 
only to the integrable cases. The equation of motion arising from the variation of 
the action ( |2.3(j| ) with respect to g gives rise to the Maxwell equation in the matrix 



potential formalism. 
Maxwell eq. : 

d{g''dg + g-'Ag) + [A, g-'dg + g-'Ag] - dA = P[T, g-'Tg]. (2.37) 
The Bloch equation again arises from the simple identity 
Block eq. : 

d{g-'fg) = [g-^fg, g-^dg + g-'Ag], (2.38) 

where we used the property [T, A] = 0. This rather abstract form of the Maxwell- 
Bloch equation will become more explicit when specific identifications of physical 
variables are made in Sec. 3. Note that the Maxwell equation is invariant under the 
if- vector gauge transformation as given in Eq. ( p.32| ), where the local function h 
now belongs to the subgroup H, while the Bloch equation is not. The integrability of 
the Maxwell equation may be demonstrated by rewriting Eq. (|2.37| ) in an equivalent 
zero curvature form in terms of the U — V pair, 

[d~U, d-V] = 0, (2.39) 



^ In a more general context, S'pot specified algebraically by a triplet of Lie groups 
F Z) G Z) H ior every symmetric space F/ G, where the Lie algebra decomposition f = g©k 
satisfies the commutation relations, 

[g , g] C g , [g , k] c k , [k , k] c g . (2.35) 

We take T and T as elements of k and define h as the simultaneous centralizer of T and 
f, i.e. h = Cg(T,r) = e g : [B , T] ^ = [B , T]} with H its associated Lie 



group. With these specifications, the action (2.30) becomes integrable and generalizes the 



sine-Gordon model according to each symmetric spaces. For compact symmetric spaces of 
type II, e.g. symmetric spaces of the form G x G/G, the model becomes equivalent to the 
type I case but with T and T belonging to the Lie algebra g. It has the coset structure 
G/H where H is the stability subgroup of T, Teg [pT], 
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and 



U = -g~'dg-g-'Ag-f3XT 
1 



V 



[2A0) 



Here, A is an arbitrary spectral parameter. This shows that the equation of motion 
becomes the integrabihty condition of the overdetermined hnear equations; 



(9- f/)^ 



{d + g-^dg + g-^Ag + PXT)^ = 
{d + A + jg-'fg)^ = 0. 



(2.41) 



The constraint equations coming from the variation of Smb with respect to A, A are 

1 



^aSmb 
SaSmb 



— I Tt{{ -dgg-' + gAg-' 
1 

2^ 



A )6A) = 



Tiii g-'dg + g-'Ag- A )5A) = 



Or, 



-Bgg-' + gAg '-A 



0, 



g-'dg + g-'Ag-A)^ = Q 



(2.42) 



(2.43) 



where the subscript h specifies the projection to the subalgebra h. It can be readily 
checked that these constraint equations, when combined with Eq. ( |2.37 ), imply the 
flatness of the connection A and A, i.e. 



[d + A, d + A]=0. 



[2.U) 



In Sec. 3, we show that various multi- level Maxwell-Bloch equations indeed arise 
from Eqs. ( |2.37|) and ( |2.43|) when appropriate choices are made for the groups G 
and H, the constant matrices T and T, and gauge fixing. 



2.3 Inhomogeneous broadening 

So far, we have obtained an action principle for the Maxwell-Bloch equation with- 
out inhomogeneous broadening. Remarkably, even in the presense of inhomogeneous 
broadening, the notion of matrix potential still persists. The inhomogeneous broad- 
ening effect, i.e. Doppler shifted atomic motions, can be incorporated beautifully via 
the U{1) vector gauge transformation. Due to the microscopic motion of atoms, each 
atom in a resonant medium responds to the macroscopic incoming light with differ- 
ent Doppler shifts of transition frequencies. Thus, microscopic variables, e.g. the 
polarization P and the population inversion D are characterized by Doppler shifts 
and they couple to the macroscopic variable E through an average over the frequency 
spectrum as given in ( p.7|) . A remarkable property of our effective field theory for- 
mulation is that it includes inhomogeneous broadening naturally only with minor 
modifications. The notion of potential variable g is again valid. In order to cope 
with microscopic motions, g becomes a function of frequency ^, i.e. g = g{z,z,C,)- 
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However, the action principle in ( |2.36D is no longer valid despite the use of the 
potential variable g. We also relax the constraint in Eq. (|2.42|) and require only 

{g-'dg + g-'Ag)^ -A = 0. (2.45) 

Then, the linear equation is given by 

L,* = [d + g-'dg + g-'Ag-^T+~XT)^ = 

where the constant A is a modified spectral parameter and becomes A + ^ in the 
absense of inhomogeneous broadening. The angular brackets denote an average over 
^ as in Eq. (|2.7|) . As in the case without inhomogeneous broadening, we make the 
same identification of the matrix g~^dg + g~^Ag — with various components 
of macroscopic electric fields which are independent of the microscopic quantity ^. 
This requires the ^-dependence of g{z,z,^) to be determined in such a way that 
g~^dg + g~^Ag — is independent of ^. It is easy to see that this requirement is 
indeed satisfied by various integrable Maxwell-Bloch systems considered in Sec. 3.0 
Note that \1/(A, z, z) is not a function of ^. The integrability of the linear equation 
(|2.46| ) becomes 



d + g-'dg + g-'Ag-iT + \T , d + (C^l 
-d{g''dg + g-'Ag - ^T) + ([T, g-'Tg]) (2.47) 







where we used the fact that g ^dg + g ^Ag — is independent of ^ and also the 
identity 

d{g-^Tg) + [g-^dg + g-^Ag, g-'fg] = 0. (2.48) 

Once again, identifying g^^Tg with components of the density matrix and Eq. ( p.48| ) 
with the Bloch equation, we obtain the Maxwell-Bloch equation with inhomogenous 
broadening. For example, we may identify E, P and D as in Eq. (p.4|) so that 
Eqs. (|2.47] ) and (|2.48| ) become the Maxwell-Bloch equation with inhomogeneous 



broadening for the nondegenerate two- level case as given in Eq. ( |2.6| ). Note that 
each frequency ^ corresponds to a specific gauge choice of the vector f/(l) subgroup. 
Therefore, in some sense inhomogenous broadening is equivalent to averaging over 
different gauge fixings of U{1) C H. This implies that inhomogenous broadening 
can not be treated by a single field theory and therefore it lacks a Lagrangian 
formulation. It is remarkable, however, that the group theoretic parametrization of 
various physical variables in terms of the potential g still survives. 



^ In the three-level system, we must take ^ — — toAi = ~tQA2. It means that in order 
to preserve the integrability in the presence of inhomogenuous broadening, two detuning 
parameters of the three-level system must be equal. 
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(a) 5f/(2)/[/(l) (b) 5f/(3)/t/(2) 



(c) SU{3)/U{2) 




(d) [SUi2)/Uil)f 




(e) SUi3)/Ui2) 




(f) 5t/(3)/[/(2) 




Jc = l 



mf, = — 1 rrif, = 1 
(g) 5f/(4)/5(f/(2) X f/(2)) 



(h) 5f/(5)/f/(4) 

Figure 1: Multi- level systems and their associated symmetric spaces 

3 Multi-level systems 

In this section, we work out in detail field theory identifications of each multi-level 
Maxwell-Bloch equations through specifying the groups G and H, the constant 
matrices T and T, and the gauge choice. Briefiy, the resulting associations with 
symmetric spaces are the following (see Fig. la-lh); 



SUi2)/Uil) 
SU{3)/U{2) 

{SU{2)/U{l)r 
SU{3)/U{2) 
SU{4)/S{U{2) X U{2)) 
SU{5)/U{4) 



nondegenerate two- level (Fig. la) 
degenerate two-level; (Fig. lb and Fig.lc) 

jb = ^ ja = 1, j6 = 1 ^ ja = 0, jfc = 1 ^ ja = 1 

degenerate two-level; j;, = 1/2 ^ = 1/2 (Fig. Id) 
nondegenerate three-level, A or V system (Fig.le and Fig. If) 
degenerate three level; ja = jc = 0, jb = 1 (Fig.lg). 
degenerate three level; ja = jc = 1, = (Fig.lh). 



(3.1) 



All of them correspond to Hermitian symmetric spaces. In the Appendix B, the 
characteristic properties of Hermitian symmetric space is used to generate infinitely 
many conserved local integrals. Our examples in Eq. (^TTJ) suggest that to each 
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Hermitian symmetric space there may exist a specific multi-level system with a 
proper adjustment of physical parameters. In particular, we could see that the multi- 
frequency generalization in a configuration of the "bouquet" type corresponds 
to the Hermitian symmetric space SU{n)/U{n — 1) for an integer n. However, for 
large n, it requires a fine tuning of many coupling constants which makes the theory 
unrealistic. 



3.1 Nondegenerate two- level system 

This is the simplest case which was originally considered by McCall and Hahn to 
explain the phenomena of self-induced transparency. It also accounts for the transi- 
tions 1/2 —s> 1/2, 1 <-> 0, 1 — * 1 and 3/2 ^1/2 for linearly polarized waves and the 
transitions 1/2 —>■ 1/2, 1 and 1 — > 1 for circularly polarized waves. We assume 
that inhomogeneous broadening is absent so that < P >= P. The Maxwell-Bloch 
equation is given by Eq. ( |2.6| ) which can be expressed in an equivalent zero cuvature 
form. 



E* -i(3\-ii) ' X\P* -D 



. (3.2) 



In order to show that this equation arises from the field theory action in Eq. ( |2.36|) , 
we take H = U{1) C SU{2) = G and T = -f = ia^ = diag(z, -i). We fix the 
vector gauge invariance by choosing 

A = i^as , A = (3.3) 

for a constant ^. Such a gauge fixing is possible due to the fiatness of A, A. Com- 
paring Eq. (|2.4CI|) with Eq. (|37^ ), we could identify E, P and D in terms of g such 



that 



g ^dg + i^g ^a-^g - i^as = { , g ^o^g = (^^ ) (3.4) 



which are consistent with the constraint equation ( |2.43| ). Note that the zero curva- 
ture equation Eq. (|2.39|) also agrees with Eq. ( |3.2|) . If we parametrize the SU{2) 
matrix g by 

g _ girders gi</3(cos6»cri-sin6lo-2)g«»7o-3 _ ( 6^**^ COS i siu (^6*^ \ ,^ 

V^sinv?e~*^ e~^*''cosv9/ ' ^ ' ' 

we recover the parametrizations of E,P and D as given in Eq. (|2.13| ) and the 
Maxwell equation becomes the complex sine-Gordon equation in Eq. ( [^.14| ). The 
potential term in Eq. ( p.36|) now changes into the population inversion D, 

^pot = / ~ 2^ = / (2-^) 
which for /3 > possesses degenerate vacua at 

{p = iPn = {n + 2)^' ^ ^ Z and 9 = 6q for 9q constant . (3.7) 

The property of degenerate vacua and the corresponding soliton solutions will be 
considered in Sec. 4. 
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3.2 Degenerate two- level system 



One of the deficiencies of the SIT model of McCall and Hahn is the absence of level 
degeneracy. Since most atomic systems possess level degeneracy, the analysis of the 
nondegenerate two-level system does not apply to a more practical system. More- 
over, level degeneracy in general breaks the integrability and does not allow exact 
soliton configurations. For example, propagation of pulses in a two-level medium 
with the transition jj, = 2 — >• ja = 2 is effectively described by the double sine- 
Gordon equation 

ddip = ci simp + C2 sin 2ip (3.8) 

which is not integrable. Nevertheless, there are a few exceptional cases which are 
completely integrable even in the presense of level degeneracy. It was shown that 
||10| , |ll| the Maxwell-Bloch equations for the transitions jb = ja = 1, jb = 
1 ^ ja = and jb = 1 — > ja = 1 (see Fig. lb and Ic) are integrable in the sense 
that they can be expressed in terms of U — V pairs. In the following, we show that 
these cases correspond to the effective theory with G = SU{3) and H = U{2) C 
G. Also, we show that the local vector gauge structure incorporates naturally the 
effects of frequency detuning and longitudinally applied magnetic field. Consider a 
monochromatic pulse propagating through a medium of degenerate two-level atoms 
in the presence of a longitudinal magnetic field. Then, the Maxwell-Bloch equation 
under SVEA is given by 



[d + l{2i + - nam)]R,m = ^ E ^'^(E Jim' " E R,,' Jl'm' 



mfi ) 



\d + Zfib(/X - /i )\R^^ = Z Y.i^'J',mRm,' - ^'*Jl'^R,n.). (3.9) 

qm 

The dimensionless quantities e'^ and R are propotional to the electric field amplitude 
E and the density matrix p, where q is the polarization index and the subscripts 
/X, /i , . . . and m,m , . . . denotes projections of the angular momentum on the quan- 
tization axis in two-level states \a) and \b) respectively.^ J denotes the Wigner's 3j 
symbols 

i;„ = (-i)— v^(_t I ;')• '"°> 

and Qai^b) is a dimensionless coupling constant of an external magnetic field. In 
general, Eq. (^.9|) is not integrable. However, with particular choices ja and jb, 



^For details of propotionality constants and their physical meanings, we refer the reader 
toRef. [§. 
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Eq. ( ^.91 ) can be recasted into the zero curvature form, or the U — V pair as in Eq. 
( |2.39| ). Specifically, for the transition jf, = 1/2 — > = 1/2 (Fig. Id), we have 



where 



U. 



± 



X 



U 



-i{x + A) 



±1* 



U+ 

u. 



i{x + A) 



V 



V+ 
V. 



{3.11] 



v+ 



R 



(ba) 
+ 2^2 



t4± 



^2^2 



In the context of field theory, we identify the U — V pair in terms of g by 

1 



U = -g-'dg - g-'Ag - f3XT, 
where the gauge choice is 

A 

and 



V 



-j9-^Tg 



(3.12) 



(3.13) 



/ 


—ix 
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—ix 





V 











ix J 



A = 0, 



T 



-T 



(T3 



(3.14) 



(3.15) 



^3, 

with the Pauli matrix era. Here, we set /5 = 1 for convenience. The resulting field 
theory is specified by the coset G/H = (Sf/(2) x SU{2))/{U{1) x f/(l)) such that 
91 



9 



92 



^3 



and 



. Note that the specific form of the identification in Eq. 



with gi,g2 C SU{2) and the two U{1) subgroups are generated by 
/O 

OJ VO (T3. 

( |3.12| ) requires gi and g2 to be SU (2) matrices as in the case of the nondegenerate 
two-level system. Thus, this case is identical to two sets of the nondegenerate two- 
level system. 



Another integrable case is for the transitions; jb = 1 ^ ja = or jb = 
In each case, the U — V pair is given by 



Ja = 1. 



u 



X 



X 



( —^i\ + i{x + y) —ie ^ 
—ie"^* |zA — ix 



—IE 





V 



1* 







i\ — iy 



2 2 

-i, y = na--i for ib = 



Ja 



2 2 

-^b - y = ^b- 2^ for jb = l ja = 0, 



(3.16) 
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and 



V = — 



2A 



V = — 



I 

2A 



Tj{ba)* 
\ -"--10 



jj{ba) 












it_l_l 


p(a) 


for 






p{a) 


p(a) 










p(H 










nil 






for 




-"■1-1 


-«!-'!- J 









(3.17) 



The gauge choice is given by 



A 




A = 0. 



(3.18) 



Thus, the field theory is specified by G/H = SU{3)/U{2) with 



T = -T 



2i 
~3 



(2 

-1 

yo 








(3.19) 



and j3 = 1. It is interesting to observe that the electric field components e'^^e ^ in 



Eq. ( p.l6|) parametrize the coset SU{3)/U{2) and the vector {e~^*,e^*)'^ transforms 
as a vector under the U{2) action. In particular, since frequency detuning amounts 
to the global [/(l)(c U{2)) action while longitudinal magnetic field amounts to the 
global U{1) X f/(l)(c f/(2)) action, the effects of both detuning and magnetic field 
to e^,e'^ can be easily obtained. 



3.3 Three level system 

The propagation of pulses in a multi-level medium with several carrier frequencies 
as given in Eq. ( p.l| ) is a more complex problem than the two-level case and in 
general the system is not exactly integrable. However, with certain restrictions on 
the parameters of the medium, it becomes integrable again and reveals much richer 
structures. Typical integrable three-level systems are either of A-type or V^-type as 
in Fig. le and Fig. If. The U — V pair for each system is essentially the same as that 
of the degenerate two- level system in Eq. (|3.16| ) . Instead of giving an explicit U — V 



pair using a density matrix, we present an equivalent expression of the Maxwell- 
Bloch equation for the A or V"-system and the U — V pair in terms of probability 
components. It is given by the the Schrodinger equation 

dci = iQiCs 

dc2 = i^2C3 

(9c3 = i{n*iCi + nic2), (3.20) 
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and the Maxwell equation 



S1C1C3 
S2C2C3, 



(3.21) 



where Si = 2tiN ^fuJi/Ti; i = 1,2 and c^; = 1,2,3 are slowly varying probability 
amplitudes for the level occupations, Qi = fiiEi/2h are the Rabi frequencies for 
the transitions i — > 3. Ei and E2 are the slowly varying electromagnetic field 
amplitudes, Hi is the dipole matrix element for the relevant transition and Ui is the 
corresponding laser frequency, and N is the density of resonant three- level atoms. If 
the oscillator strengths are equal (si = S2 = s), these equations can be put in the 
S'f/(3)/f/(2)-context with the following identifications; 



(3.22) 



and 



9 ^dg 







( * 


* * 


9 












\ * 


* * 




( 








- 
















-in* 



The gauge choice is that A = A 
matrix p, with components p^k = 



= and T = diag 
CmCl, is given by 




P 



-i9 ^Tg, T 



Finally, the U — V pair is given by 



/O 0~ 
i 
Vo 0, 



1 



(3.23) 

-i/2,i/2). The density 
(3.24) 



U = -g-'dg - sXT, V = -j9-'Tg. 



(3.25) 



This system of integrable equations exhibit many interesting exact solutions. De- 
tailed studies of this case will appear in a separate paper |2^. Recently, three-level A 
and \^-systems have received much attention in the context of quantum coherence ef- 
fects, such as lasing without inversion and electromagnetically induced transparency. 
In particular, there have been extensive studies, both analytical and numerical, on 
the propagation of matched pulses through absorbing media [^-[^. Though our 
matrix potential formulation applies only to the nonabsorbing medium case, the ex- 
act analytic solutions could provide a guideline for numerical studies in absorbing, 
non-integrable cases. It is important to note that the group symmetry persists even 



in the absorbing case which leads to interesting results pSl E3 . 



The degenerate three-level case and its integrability has been studied earlier in 



the context of the inverse scattering method 



Bloch equation formulation for the three-level case and refer the reader to Ref. 



We suppress the general Maxwell 
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for details. Here, we extend the Maxwell-Bloch equation of Ref. [IT^I to include 
a longitudinal magnetic field. Then, the Maxwell-Bloch equation describing the A 
configuration with jh = 1, ja = jc = (Fig- Ig) is given in a dimensionless form by 



de'^^ = -ipl J = 1,2, q = ±l (3.26) 



and 



id + itoihv-2A,-n,q))pl = -1(^261771^,^-8171, -sir) 

I 

I 



9 



[d - itoihv - hv - 2A2 + 2Ai))r = -I Y^ielpT - eVp\] 

g 

9^3 = -^Y.^^'jPf ~ ^TPj) 



{d + zton,{q~q))77i^^, = -^J:i^Tp'-^>T) (3-27) 

where j = 1, 2 is the amplitude of a double-frequency ultrashort pulse and 
g = ±1 denote the right (left)-handed polarization. Other variables are proportional 
to the components of the density matrix 

Pi = P-gO^^Wl-K^lX - Wit)]/Na, pI = p'Lgoexp[-i{k2X - W2t)]/Na 
ni = -pi^S/Na, 772 = -pi''S/Na, 771^^ = -p^\_jNa 



r 



-p[,o''^ exp[z(fci - k2)x - i{wi - W2)t]/Na (3.28) 



and to is a constant with the dimension of time and is the polulation density of 
the level \a). 2Ai = Wi—Wba, 2A2 = 102 — Whc measure the amount of detuning from 
the resonance frequencies. The integrability of Eq. ( p.27| ) comes from its equivalent 
zero curvature form with the 4x4 matrix U — V pair. 



U 



-Ai - i\l2x2 -iE 
-iE^ -A2 + iXl2x2 



where 



^ el e'2 J \ Pi P2 J V fii-u mil 

,'7ii r*\ . /a 0\ . / x 
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with 



a = ^ {kiv + k2V - 2Ai - 2 A2 + 4^^) 

6 = (k^v + k2V - 2Ai - 2 A2 - Ailb) 

X = !^(-3A;it. + A;2t; + 6Ai -2A2) 

y = ^-^{kiv-3k2V-2Ai + QA2). (3.30) 

Thus, in our field theory context, this corresponds to the case where G/H = 
SU{A)/S{U{2) X U{2)) and the gauge choice, 

where Ai, A2 are as in Eq. p.3(J| ) and 




(3.32) 



Similarly, we may repeat an identification for the case [^, ja = jc = 1, jfe = 0, (Fig. 
Ih) and can easily verify that it corresponds to the symmetric space SU{5)/U{A). 



4 Solitary pulses 

The lossless propagation of optical pulses in multi-level atomic media has been a 
subject of intensive study since the discovery of self-induced transparency. Most 
theoretical works on this subject have resorted to the method of inverse scatter- 
ing. The 27r-pulse of self-induced transparency and its generalizations to multi-level 
cases, e.g. simultons, are identified with "solitons" in the context of inverse scat- 
tering. Though the inverse scattering method is powerful enough to generate exact 
solutions and predict the evolution of a pulse of arbitrary shape, it does not explain 
the topological nature of solitary pulses. In the sine-Gordon limit, the 27r-pulse 
has been identified with the topological soliton of the sine-Gordon theory, which 
is stable due to the topological number conservation. The topological number is 
protected since its change costs infinite energy. The cosine potential energy in the 
sine-Gordon theory indeed measures the atomic energy through the population in- 
version. However, except for the sine-Gordon limit, such a topological treatment 
of optical pulses was not possible since field theories for more general cases were 
absent. Therefore, our field theory formulation allows a topological treatment of 
multi-level optical pulses. In the following, we show in detail that the potential 
energy term in Eq. ( |2.34D possesses infinitely many degenerate vacua and leads to 



topological solitons. In certain topological soliton is characterized by more 

than one topological number, which is a new feature of multi-level pulses. On the 
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other hand, we show that there exist also nontopological pulses which otherwise 
possess all the properties of solitons. A new, nontopological charge is introduced for 
such pulses from the "global axial [/(l)-gauge symmetry" of the field theory action 
in Eq. ( |2.36| ). Explicit nontopological soltions are constructed and identified with 
self-detuned solitary pulses. The nontopological charge measures the amount of self- 
detuning and the charge conservation law proves the stability of a nontopological 
soliton against small fluctuations. 



4.1 Potential energy and topological solitons 

The potential energy term in Eq. (|2.34|) reveals a rich structure of the vacuum 
of the theory. It is a "periodic" function in local variables. This periodicity gives 
rise to infinitely many degenerate vacua which are specified by a set of integer 
numbers. Thus, any finite energy solution should interpolate between two vacua. 
In the nondegenerate two-level case, the potential term in Eq. (|2.27|) becomes a 
periodic cosine potential in Eq. ( pl6|) and each degenerate vacuum is labeled by 
an integer n as in Eq. ( p.T] ). A soliton interpolating between two different vacua, 
labeled by and n^, as x varies from — cxo to cxd is characterized by a soliton number 
An = nh—Ua- In order to understand the vacuum structure of the potential for other 
multi-level cases, we first note that the potential term TT{gTg~^T), charaterized by 
a coset G/H, is invariant under the change g ^ gh for h & H. Consequently, we 
may express the potential term through a coset element m G G/H by Tr(mTm~^T), 
where 




cos 



- sin VB^BVB^B ^B^ cos VBW 

(4.1) 

The matrix B parametrizes the tangent space of G/H. This manifests the periodicity 
of the potential through the cosine and the sine functions. For the specific cosets, 
SU{2)/U{1), SU{3)/U{2) and SU{A)/S{U{2) x f/(2)), the relevant matrices B are 
complex- valued matrices of size 1x1, 1x2 and 2x2 respectively. Owing to the 
relation, 



B sin y/B^y/WB = sin VBB^^/bB^ B, (4.2) 
the potential term reduces to 

Tr (/ - 2 sin^ VBB^^ + Tr (j - 2 sin^ VB^b) (4.3) 

for the SU{2)/U{1) and the SU{A)/S{U{2) x U{2)) cases and 

Tr (4J - 6 sin^ V^Bt) + Tr (/ - 3 sin^ VB^b) (4.4) 

for the SU {3)/U{2) case. For a further reduction, we denote the non-zero eigenvalues 
of B'^B by 0f {i = 1, ..,T = Tank{B^B}) which are positive definite and coincide 
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with those of BB'' . In terms of 0?, the potential term takes a particularly simple 
form, 

a — 6^sin^0j, (4.5) 



where the positive constants a and b can be read directly from Eqs. (^.^l) and 



In order to check, we take for example the SU (3)/f/(2) case and choose the B matrix 
by 

5 = (-0sinr/e-*^ - cosr/e"*"). (4.6) 

Then, 



/ cos — sin sin ne '^^ 



m 



sm (J) sm rje — sin cos rje 



sin sin rje^^ cos^ r] + cos sin^ rj — cos ?7 sin r^e*^ ( 1 — cos ( 

y sin0cos?7e*" — cost] sin ?7e*"~*^(l — cos 0) sin^ r/ + cos cos^ r/ 



(4.7) 

and the potential term becomes 

TiigTg'^f) = Tr(mTm"^f ) = 6-9 sin^ (4.8) 
which agrees precisely with Eq. 



The potential term in Eq. ( |4.5| ) manifests the periodicity of the potential and 
the infinite degeneracy of the vacuum. The minima of the potential occur at 0^ = 
{rii + l/2)7r for integer n^. Therefore, the degenerate vacua are specified by a set 
integers {ni,n2, ...jUr)^ The rank r of B'^B is one for the cases of SU{2)/U{1) 
and SU{3)/U{2) and two for the case SU{A)/S{U{2) x f/(2)). Therefore, solitons 
for the SU{4:)/ S{U{2) x f/(2)) case, which interpolate between two vacua {nia,n2a) 
and (nib, 7126) with > \n2\ and rii > 0, are labeled by two soliton numbers 
Aril = nib — nia and An2 = n2b — n2a- In the following, we present an explicit 
expression for the 1-soliton carrying two soliton numbers. Consider the degenerate 
three-level system with the group structure SU{A) / S{U{2) x U{2)). For simphcity, 
we assume that the system is on resonance (Ai = A2 = 0, v = without external 
magnetic field and inhomogeneous broadening. This is equivalent to the case where 
A = A = Q'm. Eq. (|2.40|) with identifications in Eq. ( p. 301 ) in terms of a 4 x 4 matrix 



g. By applying the Backlund transformation in Ref. [Q, we obtain the 1-soliton 
solution in terms of variables as in Eq. (|4.1|), 



2 

B = -2Bq tajT^ exp{2r]z + -z + const.) = 0So (4.9) 

V 



^°In fact, in the case of multiply integer-labeled vacua, not all of them are topologically 
distinct. A similarity transformation of B^B which reshuffles the eigenvalues is a 
continuous symmetry of the vacuum, i.e. under the continuous similarity transformation, 
the potential energy does not change. For example, two vacua (711,712) and {n2,ni) are 
not topologically distinct but related by a continuous symmetry transformation. Also, 
there exists another continuous symmetry associated with the nontopological U{1) charge 
which provides an additional topological degeneracy by the identification of two vacua 
(ni, n2, n^) and (— rii, — ^2, — ri^). Thus, the topological configuration of degenerate 
vacua is characterized by {ZY/{Zr x Z2). 
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where 77 is a constant and Bq is a constant 2x2 matrix satisfying 

BoBIBq = Bo. (4.10) 
If the matrix Bq is degenerate, i.e. det-Bo = 0, it can be given in general by 

Bo = ^^('l (4.11) 



1 -|_ |q;|2 \a9i ad2 

with complex constants a and 6'i, 62 satisfying |6'ip + |^^2p = 1- The eigenvalues of 
BqBq are then zero and one. Therefore, up to a global SU (2) similarity transform of 
BqBq, this solution corresponds to the (l,0)-soliton. This solution has been known 
as a simulton in earlier literatures and its scattering behavior has been analyzed in 
detail 0. For the nondegenerate Bq, we can take Bq as an arbitrary U{2) matrix 
so that BqBq = l2x2 and the corresponding solution is the (l,l)-soliton. This is 
energetically distinct from the (l,0)-soliton and also it can not be reached to the 
(l,0)-soliton via the similarity transform since the similarity transform preserves 
eigenvalues of BqBq. Finally, physical quantities can be obtained from g through 
the identification in Eq. ( p.40|) . Explicitly, we find E, P and M in Eq. ( p.30| ) to be 

E = iBQd(j) = —2iriBosech.A 

P = -2Bq sin 2(j) = -45otanhAsechA 

M = -iV = -212x2 cos 20= -212x2(1 -2sechM) 
2 

A = 2r]z + -z + const. (4.12) 
V 

respectively. Inclusion of detuning and external magnetic effects can be done easily 
by a gauge transform; 

E ^ H^^EH2, M H^^MHi, P H^^PH2, N H^^NH2 (4.13) 

where Hi, H2 are given by Ai = H^^dHi, A2 = H2^dH2 for Ai, A2 in Eq. ( pISOl) . 



4.2 Nontopological solitons as self-detuned pulses 

Here, we address the issue of topological vs. nontopological solitons in optical sys- 
tems. In order to facilitate the problem, we first focus on the 27r pulse of the 
nondegenerate two-level system. Set /3 = 1 in Eq. ( p.6|) without loss of generality. 
Then, by using the dressing method in the Appendix A, one can obtain the 27r pulse 
solution such that 

cos 09 = , sech.(2bz — 2hCiz) 

\a - if + 



e = - tan"^ [^-^coth(2te - 2hCiz)] - 2^z + (2(a - O^i - 2^2)^ 

7] = {a-Oz+{a~OCiz, (4.14) 
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where a, b are arbitrary constants and 



where the angular bracket is as in Eq. ( p.7| ) . In terms of E as defined in Eq. ( |2.13|) , 
the 271 pulse is given by 

E = -2ib sech{2bz - 2bCiz)e-^'^'''+^^'\ (4.16) 

Note that E is explicitly independent of ^ in Eq. ([4.16|) , despite the ^-dependence 
of potential variables (f, t] and 6. This exemplifies the macroscopic nature of E as 
discussed in Sec. 2.3. In the sharp line limit of the frequency distribution f{C,) = 
— C,o), this solution retains the same form except for the change of constants Ci 
and C2, 

= {a-^oy + b^ ' = (a-^oZ+b^- ('-'^^ 



The solution in Eqs. ( [4.14[ ) and ( ^.161 ) is loosely identified with the 1-soliton in earlier 



works using the inverse scattering method. However, this does not necessarily mean 
that it is a topological 1-soliton. We emphasize that the topological distinction 
is possible only in the sharp line limit and even in that case not all solitons are 
topological solitons. For example, when a = ^, the above solution describes a 
localized pulse configuration which interpolates between two different vacua in Eq. 



(|3.71 ) such that 



Ax = -oo) = in + ^)7i, y,(x = oo) = + |^(-l)")7r. (4.18) 

Thus, it carries a topological number An = (— 1)"+^6/|6| and becomes a topological 
soliton. When a ^ C,, the solution reaches to the same vacuum as x — > ±00 since the 
peak of the localized solution does not reach to the point where cos (f = 1. That is, its 
topological number is zero. Nevertheless, it shares many important properties (e.g. 
localization, scattering behavior etc.) with the topological soliton so as to deserve 
the name, a "nontopological soliton". Note that the envelope function E, and also 
the time area of E, become complex when 07^^. But the time area of the absolute 



value of E in Eq. ( 4.16 ) is still 2tt. This suggests that we could call the solution in Eq. 



( |4.14| ) as a 2tc pulse in a broad sense, which comprises both the topological and the 



nontopological solitons as well as the inhomogeneously broadened solution. In order 
to see the physical meaning of a nontopological soliton, consider the resonant case 
where ^ = 0. In this case, Eq. ( ^.16| ) shows that the nontopological soliton (a 7^ 0) 
shifts the carrier frequency by Aw = 2a. Thus, the nontopological soliton represents 
a self-detuned 2tt pulse. It also receives a spatial modulation given by a phase factor 
exp(— 2m^/ \/a?~+l^). At the microscopic level, the maximum population inversion, 
D = cos 2ip, does not reach to 1 in the nontopological case so that the shape of is 
not of the kink-type (see Fig. 2a-2d). Note that the field intensity |£^| in Fig. 2d is 
invariantly hyperbolic secant-type independent of the value a. 
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Figure 2: Plots of (f), 6 , D and \E\ for (bright) one-soliton with 6 = 2, ^=0 
. The plot is w.r.t. x = 2bz — 2bCiZ. Real (a = 0) and dashed (a = 3) lines 
represent the topological and the nontopological solitons. 



Though a nontopological soliton can not be specified by a topological integer 
number, it carries a continuous nontopological charge. In fact, as we will show 
in the next section, the nontopological charge conservation law gives rise to the 
stability of a nontopological soliton. In Sec. 5, we show that the symmetry leading 
to the nontopological charge is "the global axial t/(l)-vector gauge symmetry" of 
the action. In the nondegenerate two-level case, this means the invariance of the 
action in Eq. ( p.l6|) under the change 



— i> ?7 + 7 for 7 constant. (4.19) 
The corresponding Noether current is given by 

COS^ 09 „ - COS^ 09 ^ , , 

J = — j^dr] , J = — ^dr] 4.20) 
sin Lf sin 

which satisfies the conservation law, 

dJ + dJ=— [cot^ ^{d + d)r]] + c— [cot^ ^dr]] = 0. (4.21) 
The corresponding conserved charges are either conserved in time, dQ'^/dt = 0, 

/+00 _ 

[cot^ ^{d + d)r]]dx, (4.22) 
-00 
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or in space, dQ^ /dx = 0, 

/ + 00 
[c cot^ If dri]dt. (4.23) 
-oo 

In the case of the nontopological sohton given in Eq. ( [4.14D , 

qS = qT = ctan-'^. (4.24) 
a -4 

The physical meaning of is clear. Consider a 2n pulse with b fixed. Since 
dr] = a — ^ = a + Wo — w expresses frequency detuning and cot^ ip is peaked around 
the soliton, Q'^ measures precisely the amount of self-detuning of a nontopological 
soliton. Stability of nontopological solitons can be proved either by using conserva- 
tion laws in terms of charge and energy as given in [Q, or by studying the behavior 
against small fluctuations. 



4.3 Stability 

The physical relevance of a topological number is that it accounts for the stability of 
solitons against "topological" (soliton number changing) fluctuations. In fact, any 
finite energy solution must approach to one of the degenerate vacua at x = ±oo 
and therefore it carries a specific topological number. Topological numbers can- 
not change during any physical process due to the infinite potential energy barrier 
between any two finite energy solutions with different topological numbers. This 
infinite energy barrier results from the infinite length of the x-axis despite the finite 
potential energy density per unit length. On the other hand, topological number 
is not useful in understanding the stability of pulses against nontopological (finite 
energy) fiuctuations. Also, the topological notion does not apply to the case with 
inhomogeneous broadening. In Sec. 2.3, we have argued that the potential variable 
g is microscopic depending on the frequency ^ while the pulse amplitude E is macro- 
scopic being independent of ^. This was apparent in the example of the 1-soliton 
solution given through Eqs. ( |4.14| )- p.l6| ). It shows that inhomogeneous broadening 
requires to be a function of "frequency ^ averaged" coefficients, while the topo- 
logical vs. nontopological nature of the soliton critically depended on ^ as in Eq. 



( 4.14 ). Thus, inhomogeneously broadened pulses do not carry topological numbers. 
In this regard, it is remarkable that the McCall and Hahn's area theorem provides 
a stability statement even in the presence of inhomogeneous broadening. In fact, 
the proof of the area theorem relied crucially on the averaging over the frequency 
^ of detuning in inhomogeneous broadening. However, one serious drawback of the 
area theorem is that it applies only to the case of real E which ignores frequency 
modulation, and it also assumes the symmetric frequency distribution. Presently, a 
more general area theorem including frequency modulation is not known. 



In this section, we attempt to generalize the area theorem to include frequency 
modulation. Though we do not have the general area theorem, we show that how 
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pulse reshaping with frequency modulation can be explained to a certain extent. 
In the nondegenerate two-level case without inhomogeneous broadening, we prove 
the stability of a 27r pulse in terms of a "modified area function" and show that 
the recovery of soliton shape is slower in the off resonant case (a — ^ 7^ 0) than in 
the resonant case. When frequency modulation is taken into account, a numerical 
work testing the pulse stability has shown that there exists a frequency-pulling 



effect [0. This effect is explained nicely in terms of the nontopological charge and 
its conservation law. Consider the 1-soliton in Eq. ( [4.14| ). Since the time area 
of complex E is not meaningful, instead we regard (f of the complex sine-Gordon 
equation as a "modified area function" . We also assume without loss of generality 
that the asymptotic time behavior of ip is given by 

(p(t = —oo,x) = — 7r/2, (p{t = oo,x) = 7i/2 for a = ^ 

ip{t = —00, x) = {p{t = 00, x) = —7t/2, for a^^. (4-25) 

Then, the modified area 

poo 

A= / 2d^dt (4.26) 



of the topological soliton (a — ^ = 0) is 27r while that of the nontopological soliton is 
zero. Consider a pulse perturbed around the 1-soliton with the boundary condition 
Lp{t = — oo,x) = — vr/2, i.e. the it is initially in the vacuum state. Near the trailing 
edge of the pulse {t » 1), the modified area function is perturbed hj 6(p = e for 
small e, i.e. (p{t » 1) = ±7r/2 + e. Then, the perturbed complex sine-Gordon 
equation for t >> 1 around the 1-soliton becomes 

The perturbation of the //-part is neglected since its contribution is of the order e^. 
This shows that if the modified area is greater than 27i (or zero) by the amount 
e > 0, then ddip < 0. Therefore, the field dip at the trailing edge tends to decrease 
along the z = x/c axis so as to recover the total modified area 27r (or zero). On the 
other hand, if e < 0, then ddip > and the field at the trailing edge increases. This 
shows that the total modified area tends to remain 2tt or zero. Moreover, Eq. ( [4.27] ) 
shows that the recovery of area is faster in the resonant case {a = C, ) than in the 
off-resonant case {a This agrees with the prediction made by a numerical work 



3T| . In fact, the recovery of the modified area can be acompanied by a stronger 
recovery of pulse shape to that of a soliton. Instead of proving this, we simply point 
out that the stability of a soliton against modified area preserving fluctuations could 
be demonstrated by modifying the Lamb's proof in terms of the Liapunov function 



32| , as well as by proving the stability of higher order conserved charges [^]. 

In order to understand the frequency modulational stability, we recall that the 
nontopological charge measures the amount of frequency self-detuning of pulses. In 
the following, we show that the stability of the nontopological charge accounts for 
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the frequency-pulling effect. From Eq. ( |4.21| ), we have 



cot' ^{d + BUztZ (4.28) 



dx 



For a 1-soliton solution, the boundary contribution is zero and is conserved in 
space. If the solution is perturbed around the soliton such that near the trailing 
edge of the pulse, 

(p{t»l,x) = ±7r/2 + e(x) 

r^{t»l,x) = («-0(t-f-^^3^^f) + 5(x) (4.29) 

for small parametric functions e{x) and S{x). To the leading order, the variation of 
then becomes 

This shows that the detuning by a higher frequency, i.e. a — ^ > reduces for 
increasing x while the lower frequency detuning does exactly the opposite. Since the 
conserved charge of the 1-soliton is ctan"^[|6|/(a — ,^)], it can be concluded that 
the absolute value of Q^{x) decreases monotonically along the x-axis. Eventually, it 
converges to a constant value of a soliton. Note that the monotonic decrease of \Q^\- 
value of a pulse is slower than that of the modified area since it is of the order e^. The 
decreasing and converging behavior of \Q^\ is in good agreement with the numerical 



work |3T[], which showed that the frequency of the optical pulse is pulled towards 
the transition frequency and reaches to a constant value along the x-axis. Thus, 
the stability of modified area and nontopological charge provides a generalization 
of the area theorem in the presense of frequency detuning in a restricted sense. A 
full-fledged generalization should include inhomogeneous broadening, in which case 
the nontopological charge conservation law breaks down. It introduces an anomaly 
term M in the current conservation law, dJ + dJ = M, for J, J as in Eq. ( |4.2CI| ) and 

M = 2 cot V5[ cos(6' - 2r]) < sm{e - 2r/) sin 2(f > 

-sm{e - 2r]) < cos{e - 2r]) sin2(^ > -(cotV^^^ + ^de)d<^]. (4.31) 



This anomaly vanishes in the sharp line limit due to the constraint in Eq. ( 2.15 ). It 
also vanishes for the 1-soliton and the charge remains conserved in this case. This 
may be compared with the conserved area of topological solitons in the presence 
of inhomogenous broadening. The area theorem of McCall and Hahn proves that 
inhomogeneous broadening changes the pulse area until it reaches to those of 2n7r 
pulses. This suggests that a generalized area theorem of pulse stability including 
frequency modulation may be proven by making use of the nontopological charge 
and the anomaly. But this has yet to be seen. 
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5 Symmetries 



One of the advantage of having a field formulation of the Maxwell-Bloch equation 
is that the field theory action reveals symmetries of the system. In this section, we 
show that our group theoretic formulation in particular reveals previously unknown 
gauge-type symmetries which have definite physical implications. Also, by using the 
group theory, we construct systematically infinitely many conserved local integrals 
of the Maxwell-Bloch equation in association with a Hermitian symmetric space 
G/H. These conservation laws can be extended to the case with inhomogeneous 
broadening without difficulty. In addition to these continuous symmetries, we show 
that the action in Eq. (|2.36 ) uncovers two types of discrete symmetries; the chiral 
and the dual symmetries. These discrete symmetries relate two different solutions. 
In particular, we show that the dual symmetry relates a "bright" soliton with a 
"dark" soliton. 



5.1 Conserved local integrals 

In the Appendix B, it is shown that the associated linear equation ( p.41[ ) in terms 
oi a.U — V pair yields exact soliton solutions through the dressing procedure. The 
same linear equation can be employed to construct infinitely many conserved local 
integrals. We first construct conserved integrals for the SU{2) /U{1) case with inho- 
mogeneous broadening and later generalize to the arbitrary G/H case. Recall that 
the linear equation for the SU{2) /U{1) case is given by 

"q^)-at)vi/ = o, (^+(x^))^ = o (5.1) 

where T = — T = ia^. We introduce the notation 







■ (Di 


Pi 




-D 



and define 



i=0 

so that the linear equation changes into 



(5.2) 



OO 1 , 

vl/exp(-AT^) = 5:-<|.,; = ) (5.3) 



E* 

and 



(^+(^* n^n-^-.-lr, $m]=0 (5.4) 



+ E = 0. (5.5) 



1=0 
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These equations can be solved iteratively in components, 

qi = ~ Esi^i) (5.6) 

Ti = -^(9r,_i + E>i_i) (5.7) 



Pi 



f Endz + 2 E / (Di-i-ipi + Pi-i-iri)dz (5.8) 

•' 1=0 •' 

i—l 

- I E*q,dz + ^ E / i-D^-l-lSl + P*_i.,qi)dz (5.9) 



1=0 ' 

together with the initial conditions: 

Po = So = -2i, ro = go = 0. (5.10) 

The consistency condition, ddpi — ddpi = 0, then leads to the infinite current con- 
servation laws, dJi + dJi = for Jj = Bpi and J, = —dpi. Or 

i-l 

idJ2{D,.i_,pi + P,_;_in) - diEn) = 0. (5.11) 

1=0 

Another consistency condition, ddsi — ddsi = 0, gives rise to the complex conjugate 
pair of Eq. (|5.11| ). A few explicit examples of conserved currents are 

Ji = -2Do 

Ji = EE* (5.12) 

J2 = 4iDi - 2PoE* 

J2 = EdE* (5.13) 

J3 = -2PodE* + 8D2 + 4:iE*Pi 

J3 = Ed'^E* + {EE*f (5.14) 

J4 = ~mD^ + 8E*P2 + AiPidE*-2Pod'^E*-2PQE*\E\^ 

J4 = Ed^E* + \E\^d\E\^ + 2E\E\^dE\ (5.15) 

Half of the above integrals have appeared earlier in As for the general G/H 
case, we introduce the abbreviation: 

S = U = g'^dg + g-^Ag-iT 

= D, + P,, (5.16) 

where in the last line, the decomposition is made according to the behavior under 
the adjoint action of T such that [T, Di] = and [T, Pi] ^ 0. Now, define matrices 
Xi and yi recursively by 

= -[T , dyi^,] - [T , S]X,.i (5.17) 
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and 

i—l 

x^ = - f sydz (A-i-i^z + p^-l-lyl)dz. (5.18) 

•' l=Q 

These matrices Xi and can be determined completely with appropriate initial 
conditions. For example, if we choose an initial condition which is consistent with 
the recursion relation for i < 0, 

Xa = I , 3^0 = 0, (5.19) 
we find for the first few explicit cases in the series, 

X^ = j 8[T , S]dz - J Dodz, y^ = -[T , S] (5.20) 

and 

^2 = J {SdS + S[T , S]X^)dz + J (-Di - DoX^ + Po[T , S])dz 

y, = -OS - [T , S]Xi. (5.21) 

These matrices give rise to infinitely many conserved local currents, 

J, = dx, = ~syi, J, = dXi = - Y,{D,_i_^Xi + Pi^i^iyi), (5.22) 

1=0 

which satisfy dJi = dJi. An explicit derivation of these currents are given in the 
Appendix B. A few examples of currents are 

Ji = £[T , £] , Ji = -Do (5.23) 

J2 = SdS + S[T , £]Xi , J2 = -Di - DoXi + Po[T , £]. (5.24) 
The first current Ji, Ji gives rise to the energy conservation law. 

5.2 Global gauge symmetries 

The action in Eq. ( |2.36| ) for the coset G/H possesses various type of gauge sym- 
metries. Since the Maxwell equation arises from the action, it also possesses gauge 
symmetries, while the Bloch equation could change under the gauge transformation. 
For example, the local if-vector gauge symmetry, as given in Eq. ( |2.32|) where the 
local function h belongs to the subgroup H, is a symmetry of the Maxwell equation, 
while the Bloch equation changes under the transformation. In fact, it was shown 
that a particular local gauge fixing accounts for the effect of frequency detuning 
and external magnetic fields. On the other hand, even after the local gauge fixing, 
there remains global gauge symmetries. For example, assume the local gauge choice 
A = A = 0. Then the action in Eq. ( p.36|) possesses the property 



SMBiLgR) = SMB{g) (5.25) 



34 



for constant matrices R and L which commute with T and T respectively. Thus, 
R is an element of the subgroup H . Under the transformation g LgR, electric 
field components the density matrix components rotate among themselves via the 
similarity transform 



g^'dg ^ R-\g-'dg)R, g-'fg ^ R-\g-^fg)R 



(5.26) 



For example, in the case of three-level A or V^-systems, the Rabi frequency Vti and 
the probability amplitude Cj are rotated by 



in which case 



h* h* 



h* 

"'12 



21 

h* 

""12 



^^2 
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"21 
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h* 
"22 







K 





h* 
"■33 




(5.27) 




(5.2^ 



Note that when i7i = and VL2 is a 27r sech pulse, the rotated Rabi frequencies are 
all propotional to the 2tt sech pulse, i.e. it becomes a simulton solution. Thus, our 
global symmetry provides a systematic way to generate simulton solutions. When 
L = R = exp(7T) G U{1), we have the global f/(l)-axial vector symmetry. The 
Noether charge of this f/(l)-invariance is precisely the nontopological conserved 
charge introduced in Sec. 4.2. Even though a general expression for the nontopo- 
logical charge should be possible, in practice it requires an explicit (noncompact) 
parametrization of the group variable g as in the case of Sec. 4.4. 



5.3 Discrete symmetries 

Besides continuous symmetries, the action in Eq. (p.36|) also reveals discrete symme- 
tries, the chiral symmetry and the dual symmetry. They are manifested most easily 
in the gauge where A = A = {). Extensions to different gauges, e.g. the off-resonant 
case which requires a different gauge fixing as in Eq. (|3.3|), can be made by the 



vector gauge transform in Eq. (|2.32| ). One peculiar property of the action in Eq. 
( |2.36| ) is its asymmetry under the change of parity through z ^ z. This is because 
the Wess-Zumino-Witten action in Eq. ( p. 26 ) is a sum of the parity even kinetic 



term and the parity odd Wess-Zumino term thereby breaking parity invariance. In 
the optics context, broken parity is due to the slowly varying enveloping approxi- 
mation which breaks the apparent parity invariance of the Maxwell-Bloch equation. 
Nevertheless, the action in Eq. ( p.36| ) is invariant under the chiral transform 



z ^ z , g ^ g (orr/^ —rj , ip ^ — V^)- (5.29) 

This may be compared with the CP invariance in particle physics. Thus, parity 
invariance is in fact not lost but appears in a different guise, namely the chiral 
invariance. This chiral symmetry relates two distinct solutions, or it generates a 
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(a) a=3, b=l 



(b) a=-0.3, b=0.1 



Cos[phi] 




Figure 3: cosy? for (a) one-soliton and (b) its chirally transformed soliton. 



new solution from a known one. For example, under the chiral transform in Eq. 
( |5.29| ), the 1-soliton solution in Eq. ( [4.14| ) in the resonant case {C, = 0) becomes 



again a soliton but with the replacement of constants a, b by 



a 



2 



62 ' a2 + 62- 



(5.30) 



This implies the change of pulse shape and the change of pulse velocity hj v ^ c — v 
(see Fig. 3). The current and the charge also change by 

J ^ - J , J ^ - J , Q^-Q. (5.31) 

It is remarkable that the velocity changes from v to c — v unlike the usual parity 
change v ^ —v. 



The other type of discrete symmetry of the action in Eq. (|2.36|) is the dual 
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13^-P , g^iag (5.32) 

where a is a constant matrix with a property, aT + Ta = 0. For example, a = ai 
of Pauli matrices in the SU{2)/U{1) case. This rather unconventional symmetry, 
also the name, stems from the ubiquitous nature of the action in Eq. ( |2.36|) , i.e. it 
also arises as a large level limit of parafermions in statistical physics and the above 



transform is an interchange between the spin and the dual spin variables [0. In 
general, the change of the sign of f3 makes the potential upside down so that the 
degenerate vacua becomes maxima of the potential and vice versa. Therefore, the 
dual transformed solutions are no longer stable solutions. This allows us to find a 
localized solution which approches to the maximum of the potential asymptotically 
(so called a "dark" soliton). In practice, the dark soliton for positive (3 can be 
obtained by replacing j3 — > —(3 , z —>■ —z in the "bright" soliton of the negative f3 
case. For example, we obtain the dark 1-soliton for the SU{2)/U{1) case as follows: 



cosipe^''^ = , :tanh(2fe + 2bCz) 



'{a-O^ + b^ ^{a-if + b^ 

e = -2{a-i){z~Cz)-2iz. (5.33) 

Fig. 4 shows profiles of a dark soliton. Note that field intensity is the same as 
that of the bright soliton in Fig. 2. However, the population inversion D for the 
dark soliton becomes inverted compared to that of the bright soliton. 
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Figure 4: Plots of (f), 6 , D and \E\ for dark one-soliton with 6 = 2, ^=0 . 
The plot is w.r.t. x = 2bz + 2bCiZ. Real (a = 0) and dashed (a = 3) lines 
represent the topological and the nontopological solitons. 



6 Discussion 



In this paper, we have introduced a potential concept to optical systems described 
by the Maxwell-Bloch equation. In terms of the matrix potential, a field theory 
action for the Maxwell equation was established where the Bloch eqaution became 
a mere identity. Various identifications of multi-level systems have been made in 
associated with specific symmetric spaces G/H and the resulting group theoretic 
properties have been used in constructing conserved integrals. The field theory ac- 
tion uncovered several new features of the Maxwell-Bloch system; gauge symmetry, 
topological and nontopological charges, self-detuning, modified area theorem etc.. 
In doing so, the introduction of a matrix potential variable g was an essential step. 
One immediate question is about the generality of such a potential variable in the 
description of nonlinear optics problems. Throughout the paper, we have confined 
ourselves only to the integrable Maxwell-Bloch equations which admit the inverse 
scattering method. Also, we have concentrated only on the classical aspect of the 
field theory which gives a semi-classical description of light-matter interaction. In 
general, the Maxwell-Bloch equation is not integrable. Even the integrable cases re- 
quire specific physical settings. For example, in the three-level system, integrability 
requires equal oscillator strengths. However, we emphasize that integrability is not 
a necessary condition for the matrix potential formalism. Even for nonintegrable 
cases, one could still solve the Bloch equation in terms of a matrix potential g, and 
express the Maxwell equation in terms of g |23|. One example is the 2 — > 2 transi- 
tion in the degenerate two-level system which is described by the double sine-Gordon 
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equation when certain restrictions are made. A more general matrix potential treat- 
ment of nonintegrable cases will be considered elsewhere. 



On the other hand, the group theoretic approach in terms of g is not restricted to 
the Maxwell-Bloch systems only. The nonlinear Schrodinger equation, which is the 
governing equation for optical soliton communication systems, can be generalized 



according to each Hermitian symmetric spaces 33]. In fact, both the Maxwell-Bloch 



and the nonlinear Schrodinger equations share the same Hamiltonian structure and 
they can be combined together. This case and its physical applications will be con- 
sidered in a separate paper. Finally, we point out that our field theory formulation 
provides a vantage point to the quantum Maxwell-Bloch system as well as the quan- 
tum optics itself. A direct quantization of the Maxwell-Bloch equation using the 
quantum inverse scattering has been made by Rupasov and a localized multipar- 
ticle state has been found and compared with a quantum soliton Our field 

theory formulation suggests an alternative, yet more systematic way of quantization 
through the use of quantum field theory. The appearance of specific coset struc- 
tures and their Hermitian properties suggests that a systematic quantization based 
on group theory is possible. Once again, this is not restricted to integrable cases and 
extensions to other quantum optical systems can be made. This work is in progress 
and will be reported elsewhere. 
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APPENDIX A: Inverse scattering method and the matrix potential 

In the following, we show that the matrix potential is intimately related to the 
dressing (inverse scattering) method and explain how to obtain exact solutions. The 
dressing method is a systematic way to obtain nontrivial solutions from a trivial one. 
In our case, we take a trivial solution by 

c/ = 1 and * = ^° = exp[-(v4 -^T+~XT)z- (^J^^ (^-1) 

Let r be a closed contour or a contour extending to infinity on the complex plane 
of the parameter A and G{\) be a matrix function on F. Consider the Riemann 
problem of \E'''G(A)(^°)~^ on F which consists of the factorization 

^°G(A)(^°)-^ = (A.2) 

where the matrix function $+(2;, A) is analytic with n simple poles /ii, /i„ inside 
F and A) analytic with n simple zeros Ai, A„ outside F. We assume that 
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none of these poles and zeros lie on the contour T and the factorization is analytically 
continued to the region where X ^ fii,Xi ; i = l,...,n. We normalize by 
= ~ Differentiating Eq. (|A.2| ) with respect to z and z, one can 

easily show that 



Since is analytic inside(outside) F, we find that the matrix functions U and 

V, defined by 

U = -9$$"^ + $(A-^r + AT)$-i - AT 

V = ~{X-^)d^^-^ + ^f^-^ (A.4) 

where $ = $+ or $_ depending on the region, become independent of A. Then, 
v[/ = satisfies the linear equation; 

{d + U + XT)^ = , {8+ /^j^^)^ = . (A.5) 

Since U, V are independent of A, we may fix A by taking A = ^. Define g by 
g = H^^^l^^^ where H is an arbitrary constant matrix which commutes with T,T 
and A. Then, U and V become 

U = g-'dg + g-'Ag-^T (A.6) 
V = g-'fg. (A.7) 

If we further require the constraint condition ( |2.45| ) on $^^|x=^ such that 

+ _ ^ = 0, (A.8) 

we obtain a nontrivial solution g and \E' from a trivial one. The nontrivial solution 
in general describes n-solitons coupled with radiation mode. If G{X) = 1 in Eq. 
( |A.2| ), we obtain exact n-soliton solutions. This formal procedure may be carried 
out explicitly for each cases of SIT in Sec. 3 and a closed form of n-soliton solutions 
can be obtained. 



APPENDIX B: Conserved local integrals 

We first review some basic facts about Hermitian symmetric space which 
are relevant for the construction of conserved integrals. A symmetric space G/H 
is a coset space with the Lie algebra commutation relations among generators of 
associated Lie algebras such that 

[h , h] C h , [h , m] C m , [m , m] C h , (B. 1) 
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where g and h are Lie algebras of G and H and m is the vector space complement 
of h in g, i.e. 

g = h©m. (B. 2) 

Hermitian symmetric space is a symmetric space equipped with a complex structure. 
In general, such a complex structure is given by the adjoint action of To on m up 
to a scaling, where To is an element belonging to the Cartan subalgebra of g whose 
stabihty subgroup is H. In our case, Tq is precisely the T-matrix given in Sec. 3. 
Namely, with a suitable normalization of T, we have 

T G h , [T , h] = and [T , [T , a]] = —a for any a G m. (B. 3) 

We decompose an algebra element G g according to Eq. ( |B. 2| ), 

i^ = ^h + ^m. (B. 4) 

Such a decomposition could be extended to a representation \E' of G = SU{n) if we 
substitute the commutator by a direct matrix multiplication and add an identity 
element ho = I to the subalgebra h, i.e. 

^ = + ^m, ^/x^/x C C '^m'^m C <ilh- (B. 5) 

In other words, any unitary n x n matrix can be expressed as a linear combination 
of SU (n) generators and the identity element ho such that 

dim h dim m 

a=0 6=1 

In order to solve the linear equation ( p. 461 ) recursively, we expand in terms of a 
power series in A, 

OO 1 

^exp(-AT^) = ET7<^'., (B. 7) 

i=0 ^ 

where 

dim h dim m 

$.= E n^ha+ E ^Lrn,^X, + y,. (B. 8) 

a=0 6=1 

With the notation in Eq. ( p.l6| ), the linear equation is given by 

{d + S)^,-[T , <f.+i] = (B. 9) 

and ^ 

+ E(A-n + p^-l-l)^l = 0. (B. 10) 

1=0 

Then, the m-component of Eq. (|B. 9| ) is 

dy,.i + sx,^i-[T , yi] = o, (B. 11) 



40 



which can be solved for by applying the adjoint action of T, 



= -[T , dy,^,] - [T , S]X,^,. (B. 12) 

Xi can be solved similarly from the h-component of Eqs. (p. 9|) and ( p. lOp such 
that 

j— 1 

x, = - f sydz (A-i-i^/ + Pi-i-iyi)dz. (b. 13) 

Finally, the conserved current follows from the consistency condition ddXi = ddXi. 

With the repetitive use of the properties of the Hermitian symmetric space, 
it can be easily checked that these conservation laws are indeed consistent with 
the equations of motion ( |2.47] ) and ( |2.48| ), which in the present convention take a 
particularly simple form: 

dS+[T , Po] = 
dDi +[£ , P^] = 
dP, + [S , A] - [T , P.+i] = 0. (B. 14) 

In general, the conserved current contains nonlocal terms. These nonlocal terms 
may be dropped out by taking the T-component of the currents. For instance, the 
T-component of the "spin-2" current conservation is 

OTiiTSdS) = dTT{TPo[T , S] - TDi) (B. 15) 

which obviously does not contain nonlocal terms. 
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